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ABSTRACT 

We describe the evolution of double degenerate binary systems, consisting of components obeying 
the zero temperature mass radius relationship for white dwarf stars, from the onset of mass transfer to 
one of several possible outcomes including merger, tidal disruption of the donor, or survival as a semi- 
detached AM CVn system. We use a combination of analytic solutions and numerical integrations 
of the standard orbit-averaged first-order evolution equations, including direct impact accretion and 
the evolution of the components due to mass exchange. We include also the effects of mass-loss 
during super-critical (super-Eddington) mass transfer and the tidal and advective exchanges of angular 
momentum between the binary components. We cover much the same ground as lMarsh et alJ (2004) 
with the additional effects of the advective or consequential angular momentum from the donor and 
its tidal coupling to the orbit which is expected to be stronger than that of the accretor. With the 
caveat that our formalism does not include an explicit treatment of common envelope phases, our 
results suggest that a larger fraction of detached double white dwarfs than what has been hitherto 
assumed, survive the onset of mass transfer, even if this mass transfer is initially unstable and rises 
to super-Eddington levels. In addition, as a consequence of the tidal coupling, systems that come 
into contact near the mass transfer instability boundary undergo a phase of oscillation cycles in their 
orbital period (and other system parameters). Unless the donor star has a finite entropy such that the 
effective mass-radius relationship deviates significantly from that of a zero temperature white dwarf, 
we expect our results to be valid. Much of the formalism developed here would also apply to other 
mass-transferring binaries, and in particular to cataclysmic variables and Algol systems. 
Subject headings: accretion, accretion disks — binaries: close — gravitational waves — stars: novae, 
cataclysmic variables — stars: white dwarfs 



1. INTRODUCTION 

Binary white dwarfs which are close enough to be driven together by angular momentum losses due to gravitational 
radiation, and perhaps additional mechanisms, will undergo a phase of mass transfer during which the ultimate fate of 
the binary is decided. There are many reasons to revisit these dynamical phases of the evolution of double degenerate 
binaries today. Chief among these is that a survey of the published literature on the subject reveals that many pieces 
of the puzzle have been explored in different degrees but we still lack a uniform theoretical understanding of the fate 
of these binaries for all possible mass ratios, orbital parameters and origins. Ideally, one would like our theoretical 
understanding to be such that given a white dwarf binary of arbitrary masses and compositions at the time that the less 
massive component gets into contact, one could reconstruct the previous evolutionary pathways and the subsequent 
evolution to merger, tidal disruption or stable mass transfer. 

The reason for this lack of uniformity is that understanding double white dwarf binaries (DWD) in this rapid phase 
of interaction can be quite complex and demanding, necessitating detailed hydrodynamics, nuclear physics, radiative 
transfer and stellar structure. Therefore it is natural that different assumptions and techniques are used when one 
attempts to answer a question relevant to different classes of phenomena. For example, the approaches adopted in 
studying the putative progenitors of supernovae of type la, or the sources of gravitational waves for LISA, or the 
progenitors of AM CVn binaries, are all very different. While we shall not attempt to cover all the rich physics that 
may be ultimately necessary to have a uniform and reliable treatment of all possible outcomes of the interaction, 
we will describe these interactions within a single semi-analytic framework, making connections where appropriate to 
well-known results already in the literature. 

Our particular motive for having developed the u nderstanding described in this paper is that our group has been 
improving and running a 3-D numerical hydrocode l)Motl et alJl2002t iD'Souza et alJ l2006) which is currently capable 
of following self-consistently the evolution of model white dwarf binaries through these rapid phases of mass transfer, 
tracking mass and angular momentum with high accuracy for over 30 orbital periods. In the course of numerous 
evolutions in which we controlled the rate of driving by angular momentum losses, or prescribed an arbitrary rate 
of pseudo-thermal expansion of the donor, it became clear t hat we needed some sim pler insight into the behavior 
of the models. This led us to extend the analytic solution of Wcbbin k fc Ibenl l)1987f) (WI) for the time-dependent 
behavior of the mass transfer by relaxing most of the assumptions made to render the problem tractable. Here, we 
retain the simplifying assumption of Roche geometry, but allow all the binary parameters to vary self-consistently. 
Thus we investigate numerically a system of first-order evolution equations for a variety of cases and recover the WI 
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solution when appropriate conditions are applied. We then extend our insight to more general cases, and comment 
on the applicability and limitations of our approach. While we are still far from the more ambitious goal described 
above, we think that the present investigation may also provide a useful framework for other workers in the field of 
binary evolution and especially for those interested in large-scale numerical simulations of these binaries. Many of the 
techniques and theoretical insights in this paper are applicable to other semi-detached binaries and may have some 
relevance to contact binaries. 

2. BASIC EQUATIONS 

Consider a binary having nearly spherical components of masses Mi and M 2 and separation a. Without loss of 
generality, we shall assume that as this binary evolves from a detached to a semi-detached state, star 2 is the one 
that fills its critical Roche lobe and becomes the donor. We define the mass ratio of the binary to be q — Mil My. 
Assuming for simplicity that the spin axes of the individual components of the binary are perpendicular to the orbital 
plane, and that the orbit is circular, the total angular momentum of the system is given by: 

Jtot = Jorb + J\ + Ji 

= M1M2 (Ga/M) 1/2 + kiMiRf u>! + k 2 M 2 Rj w 2 , (1) 

where fc, are dimcnsionless constants depending on the internal structure of the components and u>i are the angular 
spin frequencies. The first term in eq. Q represents the orbital angular momentum of the components, and the two 
other terms represent the spin angular momenta of the stars. The form of the orbital angular momentum term adopted 
above assumes the binary revolves at the Keplerian frequency Q = (GM / a 3 ) 1 / 2 , which is a good approximation if the 
stars are centrally condensed. 

For a fully synchronous configuration, it is easy to show that the total angular momentum and the total energy have 
a minimum at the same separation a m i n = [3(/i + I 2 )/ /i} 1 ^ 2 , where the Ii = k{M{R\ are the moments of inertia of the 
components, and /1 is the reduced mass. Even if the orbital frequency and the spin frequency are not synchronized 
but remain proportional to one another (u>i = /jO), there is a minimum of Jtot at some a m i n - Also the total energy of 
the system will have a minimum but in general it will not occur at a m i n .* 

For our purposes it will be sufficient to work with the approximate form of eq. given above. A more complete and 
thorough discussion of the secular and dynamical stability of polytropic binaries has been presented in two well-known 
series of papers by Lai, Rasio & Shapiro (LRS1-LRS5) and further developed with SPH simulations in papers by Rasio 
& Shapiro (RS1-RS3), who also addressed the role of mass-transfer. 

If the donor has a relatively soft equation of state and if q 1, the binary tends to become semi-detached and 
mass-transfer occurs before it falls prey to the tidal instability mentioned above. Mass transfer changes the initial 
configuration as the system evolves and may either drive the system to smaller separations and thus closer to the onset 
of the tidal instability or to larger separations and toward stability. Similarly, mass loss from the system can affect 
the dynamic stability of the system. We investigate the behavior of such systems (in particular DWD binaries) in the 
cases of conservative and non-conservative mass transfer, driven by gravitational wave radiation (GWR). 

2.1. Orbital and spin angular momentum 

A system of two p oint masses orbiting around each other, in circular orbits, radiates gravitationally 
ijLandau fc Lifshi"t3ll975(l . The loss of orbital angular momentum as a result is given by 



(£] = 32 G 3 M X M 2 M 

Ij/gwr ~ 5 c 5 a 4 ^ 
Although GWR is likely to be the only important mode of angular momentum loss from DWD binaries, for any general 
form of systemic angular momentum loss J sys , we can re-arrange eq. to obtain 

•-'orb = Jtot ~ (^1 + J2) Jotb = Jsys — (Jl + J?) i (3) 

allowing for the possibility of spin-orbit coupling of the angular momenta. The rate of change of the spin angular 
momentum of each individual star can be given as the sum of the advected or consequential angular momentum 
transport plus the effect of the tidal torques 

Ji = Miji + Ji.tid , (4) 

k = M232 + i 2 ,tid • (5) 

The term denoting the advected componen t from the donor has not been included in previous treatments concern- 
ing white dwarf donors {Marsh e t al.||20( M). but has been discussed in the case of evolved donors (see for example, 
iPratt fc Strittmatterl <|197fif l and lSavoniiei 1)19781). For this reason, a few remarks clarifying the meaning of the con- 
sequential terms are in order. In the above equations ji and j 2 indicate the specific angular momenta of the matter 

* LRS2 show that for Riemann-S and Roche-Riemann sequences dE = QdJ + AdC, where A is the angular velocity of internal motions 
and C is the equatorial circulation. Thus, as a binary evolves driven by gravitational wave radiation, circulation is conserved and the 
minima of E and J will coincide. However tidal dissipation does not preserve C and thus in general the minima will not coincide in the 
presence of tidal spin-orbit coupling. 
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arriving at the accretor and the matter leaving the donor respectively. In a conservative system these will refer to 
the specific angular momentum of the same material with respect to the center of mass of each star, but at different 
times. We assume that as the stream leaves the donor there is no back torque that could modify the spin of the donor. 
Therefore j% is entirely determined by the instantaneous conditions at the donor. However, the material traveling in 
the stream experiences a time varying torque due to the binary that changes j% at the expense of the orbital angular 
momentum alone i.e. not torquing the spins of either component. These are the a ssumptions unde rlyin g the standard 
calcul ations of j\ and the estimates of the circularization radius which go back to iFlanner jj ltl97l and lLubow fcShl 
l)1975j) . We will discuss later what are the appropriate values for j\ and j 2 , and write down just the general expressions 
here. Substituting eqs. (@J and 10) in eq. © and rearranging 



Jo 



-M 2 (ji - j 2 ) + ii.tid + J: 



2,tid 



where we have assumed conservative mass transfer. Also, we can write the tidal torque as 

kxM x Rl 



Jl.tiA — 



J: 



'Sl 

hi 



2,tid 



(fl-Ul) 

w 2 ) 



(G) 

(7) 
(8) 



where r si and r S2 are the synchronization timescales of the accretor and donor respectively (See Section 16.111 . Eq. JfJJ 
becomes 



= .1 



sys 



J orb = Jsys + M 2 (jl - h)~ 

M 2 {ji -h) 



kiMiRl 



k 2 M 2 R 2 2 



M2 

M 2 



kxMxR\ — {Sl - Ux) 

T si 



k 2 M 2 Rj^(n 

T S 2 



UJ 2j 



(9) 
(10) 



where tm 2 = ~M 2 /M 2 is the mass transfer time scale. In eq. H10[) the tidal torques have been placed inside the brackets 
although, strictly speaking, they are not consequential. In most cases encountered in cataclysmic variables (CVs) and 
low-mass X-ray binaries (LMXBs), r Sl 3> tm 2 and r S2 <C tm 2 , and thus usually lo\ 3> fl, while — u> 2 )/fl\ ~ t S2 /tm 2 - 
However, in double degenerate binaries there is considerable uncertainty about the synchronization timescales, and in 
numerical simulations of mass transfer tm 2 is orders of magnitude shorter than the typical timescales one encounters 
in long-term accreting binaries. Therefore we will retain both tidal terms an d investigate their effect on dynamical 
evolutions. Note that our eq. © is equivalent to eq. (1) of lMarsh et alJ l)2004f) . with two extra terms arising from the 
advective and tidal contributions from the donor. 



2.2. Binary separation 

We now derive the equations for the evolution of the binary separation, the radius of the donor and the Roche lobe 
radius using the above equations. From the functional form of the orbital angular momentum, we know that for a 
conservative system 



(f 



-) 

J J orb 



M 2/ 
— 1 

M 2 



q)+ 2a 



Comparing this with eq. (|5J and rearranging, we obtain 



a 

2a 



'-'sys 



M2 

M 2 



l-q- M, 



3i -32 



Jo 



•^orb^si 



^2 ) 



(11) 



The sign of the quantity in brackets determines whether mass transfer tends to expand the system and thus oppose 
the effect of angular momentum losses, or lead to enhanced contraction. The change of sign will generally occur at 
some mass ratio q a whose value we discuss below. Symbolically we may write 



a 

2a" 

(la : 



Jsys _ 
Jorb 

: 1 - M- 



Jl 



lid 



J2.tid ^2 r , 

7 W^ qa 

<Vb lvl 2 



31-32 
: J, 



orb 



(12) 
(13) 



Thus (remembering that M 2 is intrinsically negative) if q > q a , mass transfer will contribute to reducing the separation 
and, as we will see below, tend to make the binary more unstable to mass transfer. On the other hand, if q < q a , mass 
transfer will oppose the effects of driving and will tend to stabilize the binary. As q decreases on mass transfer, it is 
possible that a system which started life with q > q a may evolve to a more stable configuration, if it does not fall prey 
to the tidal instability. In eq. H12(l the tidal synchronization torques may be considered additional contributions to 
the driving and they may subtract or add angular momentum to the orbit depending on the case. Note that q a should 
be interpreted as the mass ratio at which the last term in eq. I|12|) changes sign. Its value can be estimated from 
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the initial mass ratio since it is a slowly varying function of q in the conservative case, but in general it is obtained 
self-consistently as the binary evolution is followed numerically. 

The second term in the definition of q a represents the effects of the net consequential transfer of angular momentum 
from orbit to spins. We introduce the symbol Cc = M2U1 — j2)/Joib for this term, noting that it stands for the 
consequential contribution to — dlog Jorb/dlogil^. For direct impact accretion ji — j C irc(~ b\HX) is the specific 
angular momentum carried by the stream as it hits the accretor. The approximate value in parenthesis is valid for 
a synchronous donor and & i = a(0.5 — 0.227 log q) is the distance from the center of mass of the accretor to the 
inner Lagrangian point L\ ijFrank et al.ll2002|) . With the standard definition of circularization radius ry t — R clIC /a 
(jFlannervl 119751 iLubow & Shulll975t iVerbunt fc Rappapordll98i iMarsh eTal]l200^ . M 2 ir.Wc./J n rh = f(l + q)r h ] 1/2 . 
Note, however, that the tidal coupling of the donor spin to the orbit was neglected bv IMarsh et alJ l)2004|) . and the 
consequential term proportional to j 2 ~ ^2^2 was not included either. The precise val ue of 72 depends on the deta ils 
of the flow in the vicinity of L\ in a non-synchronous donor ( See iKruszewskil l) 19631) and lCsatarvova fc Skonall ((2005)) . 
For the purposes of this investigation we will simply adopt j'2 ~ R\uj2 ■ With these definitions we may write 

Qa = 1 - Cc = 1 - [(1 + a)r h r\l - f ^ )■ (14) 

The net effect of the consequential redistribution of angular momentum in the binary depends on the sign of Cc • For 
DWD binaries, ji > 22 during the direct impact stage, and this holds even after the onset of disk accretion, when 
ji = \/GM\Ri and Cc becomes smaller but remains positive. In cataclysmic variables and low-mass X-ray binaries, 
the accretion disk returns via tides most of the angular momentum advected by the stream, the donor is almost 
synchronous, and the tidal coupling of the accretor to the orbit is very weak. However, in this case R\ <C a and thus 
it is more likely that J2 > jx- While all the additional terms in eqs. (|12l I13|) are relatively small, yielding q a sa 1, in 
some cases Cc < 0, and thus q a > 1 making these systems slightly more stable. 

The expression Ijl4|l obtained above proves very useful in the interpretation of results of large-scale numerical hy- 
drodynamic simulations of the dynamical evolu tion of binaries undergoing mass transfer, with and without driving by 
angular momentum losses (|D'Souza et al . 2006) (sec also H6.3fl . 

The Roche lobe radius for the donor is accurately given by the formula due to lEggletonl l[T? 



= Rl = 0-49g 2 / 3 
rL_ a 0.6g 2 / 3 + ln(l + q 1/3 ) > 



and so with the notation of IMarsh et alJ (|200 



where Cr L ~ 1/3 is the logarithmic derivative of tl with respect to ■ Collecting results, we get 

Rl _ 2J sys ^ ^ Jl.tid + ^2,tid _ 2M 2 ^ _ Crx _ 
i?L Jorb ^orb M2 2 

Generalizing the meaning of the symbols introduced bv IWebbink fc IbenI l)1987fl to include tidal and consequential 
terms, we can write the equivalent expressions 

R-r^ + ^w 2 (17) 

2J S ys „ Jl,tid + Jl.tid /, n \ 
ls L = —^-2— — (18) 

•^orb '-'orb 

C L = -2g„ + Cr L +2q, (19) 
where the symbol v stands for driving terms and the C for logarithmic derivatives with respect to donor mass. In the 
same spirit we write the logarithmic time derivative of the donor radius i? 2 = ^2(^2, t) 

Ro Mo 

(20) 

where V2 = (d \n R2 / dt) m 2 represents the rate of change of the donor radius due to intrinsic processes such as thermal 
relaxation and nuclear evolution, w hereas C2 usually describes changes resulting from adiabatic variations of M 2 
(jD'Antona. Mazzitelli fc RitterlH989(l . More generally, since the radial variations due to the above mentioned effects 
operate on different timescales, it is more appropriate to think of C2 as the effective mass-radius exponent, averaged 
over the characteristic timescale of mass transfer. If the donors are degenerate as in the WD case, or if thermal 
relaxation is sufficiently rapid, C2 is simply obtained from the equilibrium mass-radius relationship for the donor. But 
in non-degenerate donors, if the thermal relaxation timescale Tth becomes comparable to the mass transfer timescale 
tm 2 = M2/M2, as is thought to occur during the evolution of CVs, the effects of thermal lag in the donor radius 
become important, and C2 deviates from the equilibrium value (See Appendix IX| . 

' In the range < q < 1, the function £r L takes values between 0.32 and 0.46, and is well approximated by £r L ~ 0.30 + 0.16q for 
0.1 < q < 1 
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3. MASS TRANSFER RATE 

The following discussion of mass trans fer and its stability is rooted in a similar treatment of mass transfer under 
consequential angular momentum losses ijKing fc Kolblll99q) . For all types of donor star, the mass transfer rate is 
a strong function of the depth of contact, defined here as the amount by which the donor overflows its Roche lobe 
Ai?2 = i?2 — Rl, suitably normalized. We adopt the following expression valid for most cases of interest: 

M 2 = -M {M 1 ,M 2 , a)f(AR 2 ), (21) 

where Mq is a relatively gentle function of binary parameters, while / is a rapidly varying dimensionless function of 
Ai?2- For example, for polytropic donors with index n, f = (AR2/ R2) n+3 ^ 2 l|Paczvnski fc Sie nkiewic zlll972j) : while 
for a donor with an atmospheric pressure scale height H, f = exp [AR 2 /H) is appropriate l|Ritterlll988j) . In fact, the 
exact value of the normalization rate Mo{M\,M 2 ,d) is not important at all in steady state because the equilibrium 
rate is determined by the rate of driving: given a particular value of Mo, the depth of contact will adjust to yield the 
transfer rate sustainable by the driving. In transient situations, if the mass transfer is varying rapidly, or if the depth 
of contact becomes large, the normalization becomes relevant. 

In principle, Eqs. Q, JBJ, ljl2|l and H16(l - (|21[1 completely specify the system and can be numerically integrated, and 
we discuss some examples of such evolutions later in However, before proceeding, it is more illuminating to analyze 
the general implications of eq. (|21|l . Under the assumptions mentioned above 



dAR 2 \R2 i? L / : 

where we have neglected the slower variations of Mq with system parameters and we have also assumed that the depth 
of contact is small AR 2 « R2. Thus the mass transfer rate will be steady whenever the size of the donor varies in 
step with its Roche lobe. From Eqs. (17)-(20) and l|22|) we obtain the equilibrium mass transfer rate 

M 2 \ vl- v 2 fL - V2 



(M 2 \ 



- M 2 1 cq 2(g sta bic - q) C2 - Cl 
where 

Stable = 9a - + y, (24) 

is the critical mass ratio for stability of ma ss transfer. The alternative expression on the r.h.s of eq. (|23[) thus has 
the same form as in lWebbink fc Ibenl l)1987[) . except that here the driving and consequential terms include the effects 
of tidal coupling and direct impact accretion. Furthermore, we will allow these terms to vary self-consistently as the 
evolution proceeds (See 

For example, if we assume that the binary is synchronized, the orbit circular, the donor is a polytrope of index 
n = 3/2, and V2 = (no thermal or nuclear evolution), then (2 ~ — 1/3, Cr L ~ 1/3, and we get 



M 2 \ j sys /J< 



orb 



VM 2 /cq 2/3 -Cc- q 

which is the familiar form in the case of direct impact DWD's l|Marsh et al]l2004l) . except that here Q c is reduced by 
the contribution from the donor as given by eq. (|14(l . 

The equilibrium mass transfer rate was obtained by demanding that M 2 = in eq. I|22|) . In general, if M 2 7^ (M 2 )eq, 
one can rewrite this equation as follows: 

M 2 = -2^ J£^- (g stable - q) [m 2 - (M 2 ) cq ] . (25) 

The sign of the pre-factor on the r.h.s. is negative if q < Stable, thus M2 will tend toward the equilibrium value and 
mass transfer is stable. If q > Stable no attainable equilibrium mass transfer exists and mass transfer is unstable. 

4. ANALYTIC SOLUTIONS 

Assuming that most of the parameters characterizing the binary and the donor remain constant during (the usually 
much faster) evolution of the accretion rate, it is possibl e to obtain analyt i c solu tions for the evolution of the accretion 
rate itself. In this section we generalize the result of IWebbink fc Ibenl l|1987j) to include donors with an arbitrary 
polytropic index and with an isothermal atmosphere. These can be later compared to our integrations of the evolution 
equations in which we allow all binar y and donor parame ters to vary self-consistently and also with the results of 
large-scale hydrodynamic simulations ijD'Souza et al-lfeOOffi . 

Assuming that t he donor i n the binary can be represented by a polyt rope, the mass transfer rate is given by a 
formula derived bv lJedrzeied l)1969fl assuming laminar flow, and quoted bv lPaczvhski fc Sienkiewiczl (^972) 



(26) 
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Raising both sides of the above equation to the power 2/ (2n + 3) and differentiating, we obtain 

| ( _^ 2) 2/2„+3 = {ih) 2,2n + 3 [{v2 _ j + ^2 ( _ }] (2?) 

where we have set the factor R\,/R 2 to unity, given that in most situations AR 2 <C R 2 - The analytic solutions discussed 
here assume that the driving rate wl, the intrinsic radial variation rate v% (which includes the intrinsic thermal and 
nuclear evolution), and radial reaction exponents C2,Cl remain constant while the depth of contact changes. This is 
only approximately true; and a self-consistent solution will require numerical integrations. It is interesting first to look 
at the implications of equation (|27|) when no driving is present, because the solution is immediate and instructive. 
This is a situation we encounter in some large-scale hydrodynamic simulations of mass transfer in polytropic binaries 
ijD'Souza et al.ll2006fl . Defining a positive dimensionless mass transfer y = (—M 2 / Mo) 2 ^ 2n+3 \ and a characteristic 
time scale r = M 2 /M , eq. (|2"7Jl becomes 

d| = _C2-CL ytl+3 / 2- (28) 

The solution can be easily inverted to yield 

y(t) = y(0) [l + 2/(0)" +1 /2 (n + 1/2 )fa _ ^t/r] , (29) 

where y(0) is the initial mass transfer rate, normalized as above. This solution illustrates explicitly the role of 
(2 — Cl = 2(q s tabic — q)- In the stable case, £2 > Cl, the mass transfer decays asymptotically to zero over a characteristic 
time T c hr = r /[( n + l/2)y(0)'" +1 / 2 (C2 — Cl)], whereas in the unstable case, C2 < Cl, it will blow up in finite time equal 
to r c hr- Thus the essence of the stability of mass transfer in a binary is already contained in the simple case of no 
driving. The presence of driving exacerbates the natural instability or, in the stable case, the mass transfer settles 
asymptotically to a non-zero stable value. This is what we observe, for example, in CVs, AM CVns and LMXBs. 

Returning now to eq. (|27|l . with the same definitions as above for y and r, we obtain for the general case in which 
driving is present, 

^ = - C -^(y n+3/2 -yl i : 3/2 ), (30) 

where yeq 3 ^ 2 = — {M^) eci / Mq = {yi — ^l)t/(C 2 — Cl) is the equilibrium value normalized to M$. Note that in the 
stable case, this value is positive; while it is negative in the unstable case. Before we attempt to solve the above 
differential equation, it is clear from its form and the signs just discussed, that it describes a stable solution in which 
U ~^ J/cq when q < g s tabie- If, however, q > (/stable, the r.h.s is positive even if the mass transfer vanishes initially, and 
it just gets bigger as the mass transfer grows. Since y diverges for the no-driving case in a finite time, the driven case 
diverges even sooner. 

In o rder to obtain an analytic solution to eq. I|27|l which can be compared to the solution of iWebbink fe Ibenl 
l)1987|) . it is necessary to cast it in a slightly different form using the equilibrium rate defined in eq. Il23[l . and 
modifying appropriately the definitions of the integration variable and characteristic time. With the definitions y* = 
[M 2 /(M 2 ) q] 2/(2 ™ +3) , and 1/r* = {v 2 - ^L)(A/2/|(M 2 )eq|) 2/(2 ™ +3) , the differential equation for the evolution of mass 
transfer becomes 

^ = sgnfoO-a - sgn(y,)|^r+ 3 / 2 ) , (31) 
dr t* 

where sgn(y«) is the sign of y». Thus, for the stable case y* > 0, while < for the unstable case, and r* is defined 
positive. The general analytic solution comprising both the stable and the unstable case can be given in terms of the 
hypergeometric function, as follows 

- = y.afi(l, — \— ; 1 + - J- 7 -, 8 gn(y,)|y.|"+ 8 / a ) . (32) 
r* n + 3/2 n + 3/2 

While this solution can be easily plotted numerically, it is not possible in general to invert it to obta in y* (t) . In a few 
cases, simpler analytic forms can be obtained. For example, the case n = 3/2 yields the solution of IWebbink fc Ibenl 
(1987), and the case n = 1/2 is particularly simple: 

y» = — tan (t/r*) y* < unstable 

y* = tanh (i/r*) y* > stable (33) 



In the case of isothermal atmospheres, the mass transfer rate l|Ritterlll988|) is given by : 

M 2 = -Mo e ^- R ^' H (34) 

where, H is the scale height. This form of the mass transfer equation is much simpler to integrate than the one for 
polytropes considered above. With the same approximations and notation as in the steps leading to eq. (|30|l . and 
defining y — —M 2 /M = exp ((R 2 — R L )/H), we obtain 

ldy C2 - Cl R2 . 



y dt t H 



(y - 2/eq) ■ (35) 



1 2 3 4 5 5 1 

I L~ I 

Fig. 1. — Comparison of integrations with the analytic solution by Webbink and Iben. The mass transfer rate normalized to the initial 
equilibrium rate as a function of time in units of the initial t: Analytic (black curve) and numerical - green (q = 0.663), orange (0.613), 
cyan (0.563 same as WI), blue (0.543), magenta (0.523) and red (0.513). 



This can be easily integrated to obtain: 

y = 1 T\ ~/ \ -t/r- ( 36 ) 

1 - (i - y cq /yo)e i/T ™ 

where Tj SO = tljRiiyi — vi) is the timescale required for the driving to change the depth of contact by ~ H, and yo 
is the initial value, always positive for physically meaningful cases. In the stable case y cq > 0, and y — > y cq , while 
y C q < for the unstable case and y diverges in a finite time t^iv — Ti SO In (1 — y e q/ya)- If n o driving is present, we may 
set y cq = and integrate eq. (|35|l for an isothermal donor. The result is again simple and instructive 

y= s-r. (37) 

1 + (C2-Cl)2/o^^ 

In the stable case, for any initial mass transfer, the system will detach and mass transfer will tend to zero. In the 
unstable case, any non-zero initial mass transfer will grow and diverge in a finite time. 

5. NUMERICAL INTEGRATION RESULTS 

We now relax some of the constraints imposed in the previous section and integrate the evolution equations allowing 
the binary parameters to adjust self-consistently. Specifically, we compute the changes in the masses of the components 
(assuming conservative mass transfer), allow the binary separation to change as a result of any driving present, and 
compute the values of (2 and £l as they evolve. The values of £2 depend on the adopted mass-ra dius relationship for the 
donor . Her e we use Eggleton's interpolated zero-temperature mass-radius relationship cited bv lVerbunt fc RappaporH 
(1988) and iMarsh et alj l)2004|) . which is a good approximation for systems containing old, cold white dwarfs. It is 
possible th at the systems emerg e from a common envelope evolution with a massive hydrogen atmosphere around 
the donor ijD'Antona et al.M2006|) . or that the donor has a finite entropy such that £2 deviates significantly from the 
value obtained from the zero temperature mass-radius relationship for white dwarfs ijDelove. Bildsten fc Nelemansl 
2005). In such situations, the radial variation rate V2 is non-zero, and in fact can reduce the net driving rate vl — 1/2, 
leading to a lower value for the equilibrium mass transfer (see Eq. 1231) . Also, a mass radius exponent (£2) significantly 
different from sa —1/3 clearly affects the stability and evolution of the systems at the onset of mass transfer, and 
can lead to shrinking orbits even if the mass transfer is stable with M2 ~ 0. In our subsequent analysis, where we 
are concerned about the long term integrations of the OAE, we set z^2 = and use the zero-temperature mass-radius 
rel ationship. We note that for th e study of the onset of mas s trans fer in finite entropy systems like the ones addressed 
by iDelove. Bildsten fc Nelemansl l)2005fl fc iD'Antona et alJ l)200fil) , a more realistic model for the donor is required. 
In order to calculate £l, we need to either assume or determine from other assumptions how the mass and angular 
momentum are re-distributed in the binary during mass transfer. As we have seen in ^ this depends on the mode of 
accretion appropriate for the binary considered: does the stream impact the accretor or is an accretion disk present; is 
the mass transfer sub-Eddington and conservative, or are mass and angular momentum being ejected from the system 
following super-Eddington mass transfer. For most of the numerical integrations that follow, we use the appropriate 
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Fig. 2. — Comparison of integrations with the unstable isothermal analytic solution. The natural logarithm of the mass transfer rate 
normalized to the initial equilibrium rate for the case with q = 0.663, is shown as a function of time in units of t; so : Analytic (black curve) 
and numerical - green (q = 0.663), orange (0.613), cyan (0.563 same as WI), blue (0.543), magenta (0.523) and red (0.513). 



rate of driving by gravitational wave radiation z/l, wherever necessary we assume a constant driving rate. However, 
if one is interested in the relatively rapid phases of mass transfer that follow contact and onset, then the qualitative 
properties of our integrated evolutions do not depend strongly on these assumptions. 

5.1. Surviving unstable mass transfer 

We integrate numerically the evolutionary equations for a sample of double degenerate binaries whose mass ratios 
at the onset of mass transfer exceed the stable value. In order to mimic the conditions assumed bv lWebbink fc Ibenl 
l)1987|) in their pioneering analysis, we assume a constant rate of driving, that no mass loss from the system occurs 
even during super-Eddington phases, and that a tidally truncated accretion disk is present at all stages. For these 
choices, g s tabie ~ 0.49, and thus we expect unstable mass transfer in binaries whose initial mass ratio exceeds this value. 
In Fig. 1 we present the evolution of mass transfer for a selection of unstable binaries undergoing conservative disk 
accretion, including the specific value q = 0.563 used by Webbink & Iben (hereafter WI) to illustrate their analytic 
solution. These integrations show that the mass transfer in an initially unstable binary grows at first rapidly, peaks, 
and then evolves asymptotically toward an equilibrium rate which is also evolving as q changes. The analytic WI 
solution diverges in a finite time, while all numerical solutions reach a peak and then return to stability. The peak 
transfer decreases as the initial q approaches Stable- The equations we integrate are orbit-averaged evolution equations 
(OAE) in the sense that the rates of change of the orbital parameters are averaged over one orbital period. This 
approximation is valid as long as the evolution is not too rapid, and the eccentricity of the orbit remains negligible. 
Furthermore, our formulation does not include the full effects of tidal distortion and instability discussed by LRS. Our 
results suggest that an initially unstable binary may survive the onset of unstable mass transfer as long as the mass 
transfer does not get too big and that the separation does not get too small. See ^5.21 for a discussion of the limits of 
validity of the OAE under unstable mass transfer. 

It is also interesting to compare the analytic solution given by eq. 1)36(1 for an isothermal donor in the unstable case, 
with numerical integrations of the OAEs for various initial mass ratios. This comparison is shown in Fig. with the 
analytic unstable solution plotted as the single divergent black solid line. We elected to plot the natural logarithm 
of y, which is simply the depth of contact R2 — Rl in units of the pressure scale height H. All the integrations were 
started from the same initial depth of contact corresponding to an initial mass transfer of 10 -5 of the reference rate. 

5.2. Super-Eddington mass transfer 

Another effect may come into play as discussed by Irian fc Webbinlj l)1999j) when the mass transfer exceeds the 
critical Eddington rate and the evolution becomes non-conservative. We can incorporate this effect into our numerical 
integrations, allowing the excess mass to be blown out of the binary as a wind , carry ing away a specific angular 
momentum j w . We calculate the accreted fraction (3 following- lHan fc Webbinkl (^99), and modify the evolution 
equations as follows: 

Ji = -I3M 2J1 + Ji.tid , (38) 
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Fig. 3. — Various parameters for super-Eddington accretion with direct impact with q = 0.25, just above g s table arl d 1 = 0.35, well 
above g s tablc- The panels show from the top down: the accreted fraction f}\ the logarithm of the mass transfer rate (magenta) and the 
corresponding Eddington rate (dashed red) in Mq/ji; the mass ratio q\ and the total mass normalized to the initial value. The abscissa 
shows times in years from the time of first contact. The lower q accretes at super-Eddington rates for less time as compared to the higher 
q; and it docs so more gradually losing less mass (last panel). Even for the initially more unstable mass transfer (larger q), the fraction of 
mass lost is below 3%. 
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During the direct impact phase, it is reasonable to assume that the material blown out carries away the characteristic 
specific angular momentum of the stream j w = j±. After a disk forms, the angular momentum advected by the wind 
will depend on details of the flow beyond the scope of the present study. For example, if the material is lost mostly 
from the vicinity of the accretor, it will carry approximately the specific orbital angular momentum of the accretor. 
Since we do not know exactly how j w varies, we set j w = ji throughout. In this case eq. 1)40(1 simply reduces to eq. 
©■ While it would be possible to cast the evolution equation for the binary separation (|41|l in the same form as eq. 
((Til by defining 

(42) 



2(1 + 9) 



J, 



orb 



the explicit appearance of q above makes it obvious that q a must be calculated self-consistently during the evolution. 
Note that when (3=1, the above expression reduces to eq. ((13(1 . as it should. Fig. |3 shows two examples of evolutions 
with a super-Eddington mass-transfer phase. Because the OAE do not include tidal distortions of the components or 
dissipative effects, arising for example from friction during a common envelope phase, they always predict survival, 
no matter how high the mass transfer gets during an unstable phase. The only exception to this rule occurs if during 
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Fig. 4. — Evolution of the period normalized to the period of first contact and onset of mass transfer, for a binary with initial values 
q = 0.28, and M2 = O.125A/0. This choice yields an orbital period of 300 s at onset. The curves shown correspond to different initial 
tidal synchroniza tion times; t s1 = 500 yrs (green), 1000 yrs (red), 1500 yrs (blue) and 2000 yrs (magenta). Tidal timescales are evolved 
according to eq. 1441 . The longer the tidal timescale, there are more oscillations with higher amplitude. 



the evolution the binary separation falls below the value a m ; n at which the angular momentum and the energy of the 
binary reach a minimum. In that case, further loss of angular momentum inevitably breaks the synchronism and the 
system is secularly unstable. As the binary frequency increases, the spin frequencies of the components fall behind, 
tidal synchronization torques further reduce the orbital angular momentum, and finally a dynamical instability leads 
to a rapid merger. However, for DWD binaries the lower mass component almost always fills its Roche lobe at a 
"contact separation" a c well before this minimum is reached. Then mass transfer commences and soon the second 
term in eq. (jl 1|> rises enough to drive the binary apart saving it from a merger. If initially q > Stable mass transfer 
rises even more quickly and causes the binary to expand and recover from the instability. Clearly the OAE break 
down if during the transient the transfer is so high that a significant fraction of the donor overflows in one orbit, and 
the dissipative effects mentioned above may promote a merger even before the mass transfer gets that high. 

Suppose that during the transient, the mass transfer rises to £ times the critical Eddington rate. We can estimate 
the fraction of the donor mass transferred in a single orbit using the standard assumptions about the donor filling its 
Roche lobe, and taking for simplicity a mass-radius relationship Ri » 5 x 10 8 (MQ/Mi) 1 ' 3 cm. We find 

-M 2 P P q f M Y \~ 1/3 / M : 



. I0" a £ — — . (43) 
M 2 t \0.5M e J \0.1M Q ) v 1 

In most of the parameter space this fraction is so small that we can trust the OAEs to describe approxi mately the correct 
behav ior within the limits of the physical effects included in their derivation. Looking at Fig. 3 of IHan fc Webbinkl 
we see that £ < 10 for most cases with the exception of very rare binaries in which the accretor is already 
very close to the Chandrasekhar limit. Thus we expect a certain fraction of the binaries that come into contact with 
Q > <Zstat>ic to survive unstable mass transfer, even if it happens to be super-Eddington. However, to estimate the 
fraction that actually make it through, one would have to deal with the common envelope phase which is beyond the 
aims of the present study. 

6. APPLICATIONS 

In the above sections we have developed the basic framework for investigating the evolution of close DWD systems. 
By imposing the appropriate constraints on the OAE, we have compared our OAE integrations with analytic solutions 
illustrating the limitations of the latter. In what follows we investigate the consequences of the OAE when they 
are applied with all effects included. One immediate consequence of the OAE is the phenomenon of tidally induced 
oscillations which we describe first. Next, we apply the OAE to a grid of systems with different initial component 
masses in the Mi — M\ space and study the evolutionary outcome of each of these systems. Finally, we follow the 
evolution of a single system in order to achieve a better qualitative understanding of the results obtained for this 
system using a full 3-D hydro-code. 

6.1. Tidally induced cycles 

As a system gets into contact, the binary separation decreases at first until the mass transfer rate is high enough to 
reverse the trend in a. In unstable or near unstable cases, during this phase the mass transfer timescale tm 2 decreases 
rapidly and becomes shorter than the synchronization timescale of the accretor r Sl , allowing efficient spin-up and 
building a large asynchronism. As the separation increases, the mass transfer rate begins to fall and correspondingly 
tm 2 increases rapidly (See Fig. The synchronization timescale does not evolve as rapidly as the mass transfer rate 
and eventually tm 2 ^ T si • Now the angular momentum stored in the spin of the accretor is efficiently returned to the 
orbit. If enough asynchronism has been built up in the accretor during the spin-up phase, the additional injection of 
spin angular momentum to the orbit can cause the separation - and consequently the Roche lobe radius of the donor 
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Fig. 5. — Left Panel: Evolution of a binary with q = 0.28 with r si = 2500 yrs (Magenta) and 500 yrs (Green). Right Panel: Same for 
q = 0.24; T B1 = 2500 yrs (Magenta) and 500 yrs (Green). The time is set to when the binary gets into contact the first time. The higher 
the mass ratio, the system goes through a larger number of oscillations which are of greater amplitude. 



- to increase at a rate faster than the radius of the donor. On the other hand, the donor radius increases at a slower 
rate due to the decreasing mass transfer rate. This leads to detachment - the radius of the star cannot keep up with 
the increase in the Roche lobe radius. In general, whenever the effective driving vl — V2 > 0, the systems detach. 
Eventually, when tides synchronize the spins again, the driving will shrink the binary back into contact and accretion 
recommences, slowly at first, and accelerating as the separation a shrinks and the cycle repeats. However, as the cycles 
repeat, the deviations from the equilibrium rate decrease until the system settles to a steadily expanding behavior with 
the mass transfer following the equilibrium rate closely. In the absence of tidal effects, one would expect the orbital 
period to increase monotonically with time once equilibrium mass transfer has been established for AM CVn type 
systems. However, in the presence of appropriate tidal coupling, the systems may detach and this leads to oscillations 
in the orbital period - increasing when the system gets into sufficient contact, and decreasing when it detaches and 
the GWR dominates over the tidal terms. This behavior is evident from Figs. 0] and [31 in which we plot the orbital 
period of the binary as a function of time for different mass ratios and tidal synchronization timescales. Notice that the 
number of oscillations the system goes through is a function of both the tidal timescale and the mass ratio. For a given 
tidal timescale, the higher the mass ratio, the larger the number of oscillations a system is likely to go through. This is 
because the higher mass ratio implies that the system is more unstable which leads to a higher mass transfer rate and 
thus a higher degree of asynchronism between the accretor and the orbit. On the other hand, for a given mass ratio, 
a longer tidal timescale allows for a higher spin-up of the accretor, leading to more oscillations and higher amplitudes. 
There are, however, limits to how high or low the tidal timescale (as compared to the mass transfer timescale) can be 
in order to observe this behavior. If the timescale is too high, the spins and orbit are affectively decoupled whilst if 
the timescale is too low, the coupling is too efficient to allow for any significant asynchronism. Thus, in these extremes 
we do not observe any oscillations. 

During the evolution of the binary, the tidal synchronization times change because the masses and st ellar radii 
relativ e to the orbital separation are changing. We assume that the synchronization timescales evolve as in lCampbell 
lfT98l : 

(M 2 \ 2 (a 

and choose a normalization (t So ) that yields the desired initial timescales. As has been mentioned above, tidally induced 
detachment has implications for ultra-compact DWD systems *; in particular RX J0806 and V407. In these systems it 
is observed that t he orbital per i od is decre asing at a rate consis tent with GWR, but mass transfer is obviously underway 
in these systems l|Strohmaverl l)2002j) and lHakala et all (J2003)). Th is is at odds with the theoretical expectation that 
the orbital pe riod should increase . Build ing on ideas discussed in lLamb fc Melial (^987), and complementing cases 
considered bv iMarsh & Nelemansl (j200fJ ) . we have investigated the possibility of tidally induced detachment. We 
see that it is possible for the binary to detach, especially in the case of unstable, direct impact systems. The system 
parameters like the mass of the donor, accretor and the various angular momentum loss mechanisms are not accurately 
known for the known AM CVn systems. We have calculated the amount of time the tidal oscillations are in effect in 
the case of a DWD system after it first gets into contact by loss of GWR for different system parameters. In Tabled 
we have tabulated the relevant timescales as a function of the mass ratio q and the tidal timescale r Sl of the accretor. 
T osc = T1+T2+T3 is the timescale for which the oscillations last after initial contact. T\ represents the time when the 
binary is in contact but the separation is decreasing, T2 is the time when the system is in contact and the separation 
is increasing, whilst T3 represents the time for which the system is out-of-contact. N osc represents the number of 




i The higher the donor mass, the shorter the period at initial contact. Also, a higher donor mass makes it more likely that the system 
has an unstable mass ratio. Thus, in general, oscillations are more likely in short period systems. 
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TABLE 1 

Time spent in different regimes during the oscillation phase as a function of the mass RATIO q AND tidal timescale r ai . 
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oscillations a system encounters during its evolution. We see that for a given mass ratio, a system tends to spend an 
increasing amount of time out of contact with increasing tidal timescales. Moreover, the more unstable the mass ratio, 
the larger the number of oscillations and the timescale for which the oscillations last. 

A binary can spend a considerable amount of time in which tides are effective, especially in the case of unstable 
mass transfer. In fact, since the systems also spend quite a significant fraction of time out-of-contact, there should 
be many more systems with short periods than can be observed. However, even in the most favorable case, a given 
system spends less than 30% of its time in a regime where the system is in contact and the orbit is shrinking. Thus 
it is unlikely that tidally induced oscillations are responsible for the observations of P < for RX J0806 and V407. 
Nevertheless, the probability that we ca tch a system in contact wit h P < is enhanced significantly as compared to 
the case when there are no oscillations (Willems & Kalogcra 2005). For example, for the q = 0.26 case the system 
do es not undergo any o scillations and T\ ~ 1000 yrs for r Sl = 1000 yrs. A more promising idea is the recent proposal 
bv lD'Antona et aD l)2006|) . according to which the behavior of these systems can be understood if the donor possesses 
a substantial non-degenerate atmosphere which allows the donor to shrink as mass transfer proceeds. Under these 
conditions mass transfer tends to be more stable, at least initially, and cycles are unlikely. 

In the next section we follow a large number of evolutions with different initial conditions by integrating the OAEs for 
10 9 yr, in order to address the question of which kind of systems, and under which conditions, are likely to experience 
the cycles described here. 

6.2. Exploring Evolutionary Outcomes 

We explore the parameter space for DWD binaries by investigating the different types of binary evolution that 
can occur under a variety of initial assumptions. Our procedure consists of selecting the initial masses for binary 
components and evolving them under driving by gravitational radiation for 10 9 yr. In the most general case we include 
tidal coupling between the orbit and both components, we allow mass loss in super-Eddington cases, and we include 
the transition from direct impact to dis k accretion. We also choose the tidal synchronization timescales for either the 
accretor alone, as in lMarsh et alJ l)2004[) or for both components allowing even the donor to become non-synchronous. 
As a check, we have evolved a grid of models suppressing the tidal and advective terms from the donor and choosing 
the same synchronization timescale as lMarsh et alJ (|2£)04). Under these assumptions our results are indistinguishable 
from theirs. 

Before presenting the results of the evolutions, we investigate the expected behavior of the systems, focussing on 
the stability limits and whether the mass transfer is super-Eddington or not. In Fig. 0we plot the stability limits 
for two cases: one with the donor spin properly accounted for (blue line) and the other with the donor effects ignored 
(magenta line). In addition for each case, we plot the locus of points for which M2 eq « MEdd> i.e., the locus of 
points that defines a transition from super-Eddington accretion to sub-Eddington accretion. Note that these curves 
represent the equilibrium mass transfer and Eddington rates of the system at initial contact where we assume that 
the systems are synchronous and thus the tidal terms are zero. For direct impact systems the accretion rate is always 
super-Eddington if Mi > 0.21M Q (0.17 M Q if donor terms are neglected). This is because of two reasons - i) in 
general, the higher the donor mass, the higher the mass ratio and thus the systems are closer to instability, which in 
turn implies a higher M2 eq and ii) the higher the accretor mass, the lower the threshold for super-Eddington accretion 
IHan & Webbinkl li999) . However, as can be seen from Fig. the transition of the systems from direct impact to disk 
leads to changes in the stability properties of these systems - they tend to be slightly more stable because the loss 
of orbital angular momentum to the accretor spin in disk systems is smaller than in the direct impact case. This is 
reflected in the slight upturn in the stability curve around Mi ~ 1 M Q (0.85 M Q if donor terms are neglected). The 
result of this is that the equilibrium mass transfer rate for the disk systems is lower than it would have been for that 
same system if it were a direct impact system. Consequently these systems tend to undergo sub-Eddington accretion, 
and the locus of the transition between super and sub-Eddington accretion systems follows the stability curve for both 
cases. However the disk transition 'saves' only a few systems because DWD binaries with M-i > 0.25M Q , (0.23 Mq if 
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Fig. 6. — Mass-transfer stability limits (dash-dot lines) and super (above) and sub-Eddington (below) accretion boundaries (solid lines) 
with (bl ue) and without (m agenta) the donor terms included. The thin dashed black line divides direct impact systems from disk accretion 
systems (Marsh & Stccghs 2002). Because the transition from direct impact to disk accretion makes mass transfer more stable, both the 
stability limits and Eddington accretion rate boundaries follow closely the locus of that transition toward higher donor masses (see text for 
details). 



donor terms are neglected) are all super-Eddington initially. 

We present now the results of the numerical integrations of the OAE for a grid of models in the M2 — M\ parameter 
space. In Fig. [7J we plot the results for two extreme values of the accretor tidal synchronization timescales - one 
with an extremely large tidal timescale (10 15 yrs) corresponding to inefficient coupling and one for a short timescale 
of 10 yrs corresponding to highly efficient spin-orbit coupling. Firstly, we note a significantly increased region (as 
compared to when the donor terms are ignored) of parameter space corresponding to sub-Eddington mass transfer 
and probable survival as a long-term, stable mass transfer binary of the AM CVn type. Also, from our discussion 
above, we expect that the case with inefficient tidal coupling (left panel in Fig. [7J to match the curves in Fig. The 
transition from super to sub-Eddington accretion overlaps the stability boundary until Mi « 0.21 Mq after which 
it follows the stability curve defined by the direct impact to disk transition, in accordance with our expectations. 
On the other hand, since the tides almost always act to stabilize the system by effectively reducing the driving rate, 
a simplistic analysis based on comparing the initial equilibrium mass transfer and Eddington rates has only partial 
validity. While the threshold for super-Eddington accretion remains unaffected by the tidal coupling, the mass transfer 
rate is significantly lowered due to the reduced driving rate, and can be below the Eddington accretion rate. Systems 
which have super-Eddington accretion rates when the tidal coupling is inefficient, can thus accrete at sub-Eddington 
rates if the tidal coupling is strong. As a consequence we expect the domain of sub-Eddington accretion to extend to 
higher donor masses, and on comparing the left and right panels of Fig. 0this is what we observe. Finally, the locus 
of systems that undergo oscillations in their separation and binary period is illustrated in Fig. |H1 Whether a system 
undergoes oscillations or not, depends primarily on two factors: the accretor tidal synchronization timescale r Bl and 
the mass ratio q. As can be seen from Fig. |SJ oscillations are seen to occur only in a certain domain around the 
stability curve, and that domain decreases with increasing tidal synchronization times. This is due to the fact that, 
for higher donor masses, the mass transfer timescale tm 2 is quite short and if the tidal timescale r si is much longer 
than tm 2 , the spin and the orbit are effectively decoupled. As a result there is minimal return of the spin angular 
momentum to the orbit and consequently, there are no oscillations for these long timescales. In fact, referring to Fig. 
|H1 we observe that systems with high donor mass (M2 > O.35M0) do not undergo any oscillations at all. Again, this 
is because the mass transfer rates are high in this domain and thus the tidal timescales (considered here) are much 
longer than the mass transfer timescales. Consequently, the radius of the donor keeps up with the increase in the 
Roche lobe radius throughout the evolution, and the systems stay in contact. 

On the other hand, systems to the bottom right (the ones with low donor mass and high accretor mass) are stable 
systems and are more likely to be disc-systems. Thus the accretor is not spun up as much as in the case of less stable or 
unstable systems. Moreover in the disc systems, tides are extremely efficient in returning angular momentum back to 
the orbit. Both these factors conspire to decrease the level of asynchronism achieved by the accretor and consequently, 
these systems do not undergo oscillations. 

We have also studied the effect of the donor's tidal synchronization timescale on the domain over which systems 
undergo tidally induced oscillation cycles. In most cases, the level of the donor's asynchronism is rather low. Thus the 
magnitude of the term associated with the tidal effects of the donor in the driving rate vl (eq. (|19j) 'l is much smaller 
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Fig. 7. — Evolution for 10 9 yrs for an initial t s1 of 10 15 yrs (left panel) and for t s1 of 10 yrs (right panel). This synchronization time 
is for the initial configuration and evolves according to eq. 1441 . r S2 is calculated based on whatever t so is required to obtain the desired 
value of t Si . The red symbols represent super-Eddington accretion, the green are sub-Eddington, the pluses (+) and hollow diamonds 
(0) represent systems with total mass below and above the Chandrasekhar limit, respectively. Among the latter, asterisks over diamonds 
indicate systems in which the accretor does not reach the Chandrasekhar limit in 10 9 yr. The blue dash-dot line is the initial stability 
boundary with all donor effects included. The magenta line is the stability boundary without these effects, shown here for comparison. 
Note the transition in the super and sub-Eddington accretion rate around an M2 of 0.2 Mq to the latter stability curve. 



than the corresponding term for the accretor. As a result, the donor synchronization timescale has a limited impact 
on the domain over which the systems oscillate. 
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Fig. 8. — Systems which undergo 'oscillations' at least once during their evolution for r 81 = 500 yrs (pluses, 4), 2500 yrs (crosses, x) and 
5000 yrs (open circles, o). Red indicates super-Eddington transfer during any part of the evolution, whereas green indicates sub-Eddington 
transfer throughout the 10 9 yr evolution. The magenta and blue lines are the stability limits as in previous figures. r S2 is held at a constant 
100 yrs. 



6.3. Comparison with hydrodynamic simulations 

One of the stated objectives of this paper is to develop a theoretical framework to interpret and analyze results o f 
large-scale, self-consistent, 3-D hydrodynamic simulations of binaries undergoing mass transfer Ip'Souza et al.ll2006|) . 
To this end, we apply the same initial conditions to our OAE as have been to the various runs carried out by D'Souza 
et al. We have the tidal normalization factor (r So , see eq. (03J) and the mass transfer rate scaling (m) - such that 
M2 = — mMo/(A) - as 'free parameters', which we adjust so as to obtain as close a match to the hydro-runs as we can. 
In the particular run shown in Fig. t Sq = 0.75 and m = 35.0. This choice is not unique - indeed, we get reasonable 
'fits' even for other c ombinations of these 'free parameters'. It should be noted that in the numerical simulations of 
iD'Souza et all (2006) the minimum resolvable mass transfer was on the order of ~ 10 5 M/P, which translates for 
short period AM CVn binary parameters to ~ 10 4 A/Edd! Also, in the case of the hydro-runs one observes severe 
distortion of the accretor and the formation of an accretion belt around the accretor towards the end of the evolution. 
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Fig. 9. — Comparison with some of the numerical simulations (q = 0.5 run) in D'Souza ct al. 1 200(f). Three simulations were performed 
for the same initial conditions but the binary was driven by angular momentum losses at the level of 1% per orbit for different times in 
order to achieve increasing depths of contact: Q0.5-a (blue; driven initially for 2.7 orbits), Q0.5-b (green; driven for the first 5.3 orbits) 
and Q0.5-C (red; driven throughout). The solid black curves show the accretion rate in donor masses per period, the binary mass ratio, 
the separation normalized to the initial separation, and the spin angular momenta as predicted by the OAEs, while the dotted color curves 
show the same quantities as derived from the results of the simulations. We arbitrarily change t b; of the donor and accretor to match the 
Q0.5-a run and then predict the outcomes of the other simulations. Here, r si ~ 150P and t S2 ~ 3.5P initially. 



These features cannot be easily incorporated in the OAEs, and so the later stages of evolution especially in the case 
of Q0.5-C, cannot be properly represented in the OAE. An added complication is that for the hydro-runs, the driving 
was cut off after the systems were thought to have reached a deep enough contact. Since the effective density levels, 
especially near the edge of the stars in the 3-D numerical model differs from an ideal n = 3/2 polytrope due to the 
finite numerical resolution of the code, the depth of contact achieved after a certain amount of driving for a certain 
period of time is not necessarily the same for the hydro-runs and the OAE runs. This is especially true for Q0.5-a, 
because it is the one which is most sensitive to the depth of contact at the instant the driving is cut-off. For runs 
Q0.5-b and Q0.5-C, the depth achieved is deep enough to make small differences between the hydro-runs and the OAE 
unimportant. We are working on another set of runs in which the driving is not cut off and the systems are driven at 
slightly lower (and more realistic) rates. We hope that this will eliminate another source of discrepancy between the 
hydro-runs and the OAE. 

Despite the above mentioned shortcomings, the OAEs do reproduce reasonably well the behavior of the binaries 
we have studied. One notices that the hydro-runs have a gentler slope initially which progressively gets steeper as 
compared to the OAE runs. This, we believe, is a consequence of the complicated fluid flow around the L\ point and the 
distortion of the donor star. Moreover, during the initial stages of the evolutions, the 3-D hydrodynamic simulations 
are subject to numerical noise which is not the case for our numerical integrations. The relative significance of this 
noise diminishes as the mass transfer rate inc reases during the evolution. Also, the epicyclic motion that one encounters 
in the hydro-runs (see lD'Souza et all l|200fifl for details) cannot be reproduced in our results, since we assume circular 
orbits. Thus the behavior of the OAE is much smoother than the numerical hydro-runs with no abrupt changes in 
slopes of the various parameters. 

We conclude from the above discussion and Fig. that the OAEs confirm that: a) tidal effects play an important role 
in the numerical simulations of the binary, b) direct impact accretion is an important effect and can lead to significant 
spin-up of the accretor at the expense of orbital angular momentum, and c) the OAE prediction that systems that 
are initially unstable can indeed survive mass transfer seems to bear out in the hydro-runs despite the rather extreme 
conditions to which the binaries in the hydro-runs are subjected. 

In addition to the Q0.5 runs presented above, iD'Souza et alJ 1)200(111 also present other runs which result in the 
disruption of the binary. As can be seen from their Fig. 3, 6 and 10, the donor star becomes increasingly distorted 
towards the end of these simulations. The OAE, at least in their present form, are not capable of accounting for 
the distortion of the components and the consequence of these distortions on the fate of the binary. Additional 
hydro-runs with different values of the mass ratio q and low er driving rates v (though th ese rates are still orders of 
magnitude higher than realistic values) are being carried out ijMotl. Tohline fc~Fr ank 2006). In principle this can help, 
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for example, in determining a limit on the mass ratio q (> (/stable) above which the tidal disruption of the binary is 
likely. However, other physical effects not yet included in the 3-D simulations, such as radiation and the formation of 
a common envelope, are more likely to determine the fate of such systems. 

The tidal time scales that most closely account for the behavior observed in the simulations will serve as a measure 
of the numerical dissipation present in the simulations. We will present elsewhere a more comprehensive study of the 
numerical dissipation and its dependence on spatial resolution. 

7. DISCUSSION 

We have re-examined the standard circular orbit-averaged equations (OAEs) that describe the evolution of mass 
transferring binaries allowing for advective and tidal exchange of angular momentum between the components. We 
found that the mass transfer stream issuing through the L\ point has two effects in the internal redistribution of 
angular momentum in the binary: 1) it delivers a specific angular momentum j c - 1IC to the accretor spinning it up at 
the expense of the orbital angular momentum; and 2) it reduces the spin angular angular momentum of the donor 
by a specific amount ~ i?2^2, and ultimately couples to the orbit via tides. In the cases examined in this paper, the 
effect of this additional term is mildly stabilizing. For example, with parameters thought to be appropriate for the two 
short-period, direct impact binaries V407 Vul and RX J0806, the additional donor term is ~ 0.1 (see eq. JSJ), so that 
th e net effect of the c onsequential terms (accretor plus donor) ss —0.3, is still de-stabilizing but less than estimated 
bv lMarsh et alJ l)2004j) . and q a w 0.7. Its effect on the evolution of other systems remains to be explored further, but 
in systems driven by magnetic braking, it is likely to be relatively minor since it is relatively smaller and would be 
masked by the magnetic torques acting directly on the donor. 

We have extended analytically and numerically our understanding of the evolution of stable and unstable mass 
transfer in semi-detached binaries, with sp ecial emphasis on DWD binaries. In particular we have extended the 
analytic solutions of the type discussed bv iWebbink fc Ibenl l)1987[) to other polytropic indexes and the isothermal 
case. The analytic solutions predict that if q > (/stable at contact, the mass transfer rate diverges in a finite time 
implying the catastrophic merger of the two components. The OAEs on the other hand, always predict that a binary 
undergoing unstable mass transfer would survive after a phase of rapid mass transfer which in many cases would reach 
super-Eddington levels. This is because, unlike the analytic case, in the numerical integration of the OAE we allow the 
binary parameters to evolve self-consistently throughout the evolution. Thus, even if initially q > (/stable, at some point 
in the evolution q < (/stable, and the systems settle to the equilibrium mass transfer rate corresponding to the current 
values of the system param eters. While we do inco rporate mass loss due to super-Eddington accretion, following the 
treatment along the lines of lHan fc Webbinkl l)1999l) . clearly the OAEs must break down if a common envelope forms. 
The treatment of common envelope evolution is beyond the scope of the present paper, but we suspect that many 
systems previously considered doomed to a merger would actually survive provided that the mass transfer peak is not 
too high. This point clearly needs further investigation and probably requires 3-D hydrodynamic simulations with the 
inclusion of radiative effects. 

One interesting consequence of the tidal coupling of the accretor to the orbit is the appearance of mass transfer 
oscillation cycles occurring for a limited time after the onset of mass transfer. The likelihood of these cycles is higher 
in situations where the mass transfer is high, or rises to a high value rapidly. Therefore they tend to occur near and 
around the stability boundaries. Note that all the systems that undergo oscillations are direct impact systems, at least 
for the system parameters and synchronization timescales we have considered. Thus it is highly unlikely that systems 
that have accretion disks will undergo the tidally induced oscillation cycles in the case of DWD binaries. 

The presence of a m assive (~ 0.01M Q ) non-degenerate atmosphere on the donor at the onset of mass transfer 
(p'Antona et alJ l20061 can affect the stability and evolution of these systems. For example, a large and positive (2 
would imply (/stable > 1, and has significant implications for cycles, stability boundaries and super-Eddington mass 
transfer. The full consequences of these circumstances remain to be explored further. 

DWD systems are some of the most common compact systems in the galaxy and are of particular importance for the 
space based gravitational wave detector LISA. AM CVn systems are guaranteed sources for LISA and the knowledge 
of possible evolutionary trajectories is valuable. The framework we have outlined in this paper can be used to generate 
templates for sho rt period DWD's in general an d A M CVn system s in pa rticular. Similar work has been done already 
(see for example, iKopparaou fc Tohlind l)2006|) and IStroeer et ahl (|2005f) ^: but the effects of the tidal couplings and 
the advective term associated with the donor spin remain to be incorporated into future studies. 

The LSU theory grou p has performed a number of large-scal e 3-D numerical simulations of interacting binaries with 
polytropic components ijMotl et al.ll2002t iD'Souza et alJ l2*006) . These simulations did not include enough physics to 
tackle the common envelope evolution, but they should be viewed as the first steps toward that goal. In the meantime, 
we have used the OAEs with s uitably adjusted tidal coupling time scales to analyze and interpret the results of some 
of the simulations described by ID'Souza et al.l l)2006j) . The mass transfer rates that these simulations can resolve are 
much higher than the Eddington critical rate and probably much higher than the rates likely to arise during the onset 
of mass transfer in most realistic cases. Nevertheless they describe correctly the dynamical aspects of the mass transfer 
and tidal interactions under these conditions. Comparing the predictions of the OAEs with the simulations, we were 
pleased to find that they predict the outcome of the simulations reasonably well. 
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Joel E. Tohline, Patrick M. Motl and Ravi Kopparapu for helpful discussions. Finally, we would like to acknowledge 
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APPENDIX 

THE EFFECTIVE MASS-RADIUS EXPONENT C2 

We derive here a simple analytic approximation to the effective mass-radius exponent when the response of the donor 
is a combination of the adiabatic an d thermal adjustments to mass loss. Our starting point is the same as in §3 of 
ID'Antona. Mazzitelli k, Ritterl ljl989fl . namely 

dlni? 2 fdlnR 2 \ f d\nR 2 \ M 2 

+ Qs-r-r, (Al) 



d* V dt J th V dt J nuc -M 2 

where the first two terms on the r.h.s. represent the effects of thermal relaxation and nuclear evolution, and Cs is 
the purely adiabatic mass-radius exponent. The thermal relaxation term may arise as a result of initial conditions, 
nuclear evolution and mass transfer, and it is usually not possible to disentangle these effects if they operate on similar 
timescales. However, an approximate description of the radial evolution can be obtained by viewing it as the result of 
the superposition of thermal relaxation of initial conditions and nuclear evolution, plus a thermal adjustment to mass 
loss, as follows 

dlni? 2 , j. Mi , A „s 

-5- = "* + 4 + ( ) 

where v 2 is the superposition of intrinsic thermal and nuclear evolution, while v' 2 stands for the rate of thermal radial 
reaction to mass loss. Our goal is to combine the last two terms on the r.h.s. by absorbing the thermal adjustment to 
mass loss into an effective mass-radius exponent. We write 

C2 mt^ +Cs m2 (A3) 

where (2 is the effective mass-radius exponent we seek. With the above interpretation, the term v 2 in eq. I|2U|I 
represents intrinsic thermal and nuclear processes which may operate on timescales that differ from the standard 
thermal relaxation rate. This approach is consistent only if v 2 « C, 2 (M 2 /M 2 ). Therefore, in what follows, we shall 
only consider thermal relaxation in response to mass loss, and set v 2 = 0. As a consequence of the mass loss, the 
donor's radius R 2 will differ from the equilibrium radius corresponding to its instantaneous mass R cq (M 2 ). With these 
definitions we write 

R 2 _ R cq (M 2 ) -R 2 M 2 
R 2 R 2 r M 2 

where we have adopted a simple linear approximation for the thermal reaction term and t' is the corresponding 
time scal e. The secular evolution of the binary takes place on a mass-transfer timescale tm 2 = —M 2 /M 2 . Differentiating 
eq. (IA4|) with respect to time, we get 

a? i„ e>_ 1 z? ft t. \ r 

(A5) 

If the effective mass-radius exponent is to have any meaning, it must not change much over the evolutionary phase 
one is considering. Thus, we require 
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= -b/rt . (A6) 



dt 2 

Finally, setting R cq = R 2 in eq. QA5J1 . and solving for C, 2 , we obtain 

> Coq + (s T> /l~M 2 I a -s 

6 = l + rVr M2 • (A7) 

This expression shows that if the evolution is much slower than the thermal relaxation (r' <C tm 2 ) > the donor radius 
follows the equilibrium radius closely, whereas if mass transfer occurs rapidly, the donor reacts adiabatically. The above 
discussion may be applied to cataclysmic variables to describe approximately how the donor becomes increasingly 
oversized as the orbital period decreases because v 2 w 0. If relaxation from initial conditions or significant nuclear 
evolution is taking place on timescales comparable to either r' or tm 2 , the above discussion is strictly not valid. 
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